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agreement  between  the  Navy  and  the  NSF,  whereby  two  3\ACPI  probes  are  being  developed  in  parallel  under  joint 

funding  from  both  agencies.  Both  the  Navy  and  NSF  agreed  to  this  arrangement  to  facilitate  development  of  the  3V-CPI 
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1.  Introduction 


SPEC  Incorporated  was  contracted  by  the  Office  of  Naval  Research  to 
develop  innovative  new  technology  to  measure  the  size,  shape  and 
concentration  of  water  drops  and  ice  particles  in  clouds.  This  effort  included  the 
main  SBIR  Phase  II  contract  and  two  options.  The  principal  objectives  of  the 
main  SBIR  phase  II  contract  and  its  two  options  are: 

Main  Contract  Development  of  a  2D-S  (Stereo)  cloud  particle  imaging  probe 

Option  1:  Determine  the  feasibility  of  designing  and  developing  a  3V-CPI 
probe. 

Option  2:  Design  and  fabricate  a  3V-CPI  probe  for  use  on  research  aircraft. 

One  result  of  the  Phase  II  research  was  the  development  of  a  new  cloud 
particle-imaging  probe,  the  2D-S  (Stereo)  probe,  which  has  been  installed  on  four 
research  aircraft  and  flown  in  four  major  field  campaigns  sponsored  by  the 
National  Science  Foundation  (NSF)  and  the  National  Aeronautics  and  Space 
Administration  (NASA).  A  description  of  the  2D-S  probe  and  scientific  results 
from  field  campaigns  are  the  subject  of  two  refereed  journal  publications  (Lawson 
et  al.  2006;  Lawson  et  al.  2007). 

A  second  result  is  the  ongoing  development  of  another  innovative  cloud 
particle-imaging  probe,  called  a  3V-CPI,  which  is  essentially  a  combination  of  a 
2D-S  and  the  SPEC  cloud  particle  imager  (CPI).  The  3V-CPI  is  being  developed 
under  a  joint  agreement  between  the  Navy  and  the  NSF,  whereby  two  3V-CPI 
probes  are  being  developed  in  parallel  under  joint  funding  from  both  agencies. 
Both  the  Navy  and  NSF  agreed  to  this  arrangement  to  facilitate  development  of 
the  3V-CPI  under  available  funding.  The  3V-CPI  is  a  very  complex  instalment 
that  is  being  developed  for  a  new  NSF  research  aircraft,  a  Gutfstream  V  (G-V), 
designated  HIAPER  (High  Altitude  Platform  for  Environmental  Research). 
Subsequent  to  award  of  Phase  II  Option  2,  the  NSF  decided  to  change  their 
approach  to  modifying  the  G-V,  and  instead  of  seeking  public  category  or 
restricted  category  designations,  the  NSF  sought  to  have  the  research 
instrumentation  certified  in  transport  category.  The  distinction  is  that  transport 
category  is  the  same  as  the  normal  category  for  transporting  passengers,  while 
public  and  restricted  category  limits  occupants  to  only  crew  members.  The  FAA 
requirements  to  have  the  research  instrumentation  certified  in  transport  category 
are  much  stricter  than  public  or  restricted  categories  and  require  more  extensive 
design,  tests  and  documentation.  The  decision  to  certify  the  G-V  in  transport 
category  significantly  lengthened  the  timeline  required  to  design  and  fabricate  the 
3V-CPI.  Thus,  the  3V-CPI  was  not  completed  at  the  (mandatory  5-yr)  expiration 
date  of  this  SBIR  Phase  II  contract  and  will  be  retained  by  SPEC  after  the  close 
of  the  contract  for  final  assembly  and  testing. 
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2.  The  2D-S  (Stereo)  Probe. 

The  2D-S  (stereo)  probe  was  specifically  designed  to  overcome  limitations  of 
existing  probes  that  image  cloud  particles  in  the  size  range  from  30  pm  to  about 
1  mm.  Here  we  discuss  briefly  the  limitations  of  current  imaging  probe 
technology  and  how  the  new  probe  overcomes  these  limitations.  Errors  due  to 
electro-optical  limitations  that  are  inherent  in  the  design  of  the  Particle  Measuring 
Systems  (PMS)  model  2D-C  and  Droplet  Measurement  Technologies  (DMT) 
model  CIP  imaging  probes  have  been  well  documented  in  the  literature.  Korolev 
et  al.  (1991,  1998)  and  Strapp  et  al.  (2001)  have  shown  that  large  counting  and 
sizing  errors  occur  with  the  PMS  2D-C  particle-imaging  probe  when  drops  pass 
outside  the  optical  depth  of  field  (DOF).  For  example,  a  100  pm  particle  can  be 
oversized  by  nearly  a  factor  of  two,  undersized  by  a  factor  of  four  or  missed  all 
together.  Also,  small  (e.g.,  25  pm)  particles  have  a  very  small  DOF  (about  1  mm 
or  1/60th  of  the  distance  between  the  probe  arms)  and  very  large  uncertainties  in 
sample  volume. 

Additional  counting  and  sizing  errors  occur  due  to  the  relatively  slow  time 
response  of  the  PMS  2D-C  and  DMT  CIP  probes,  resulting  from  the  finite  time 
response  of  the  photodiode  array  and  first  stage  amplifier.  Baumgardner  and 
Korolev  (1997),  and  more  recently,  Strapp  et  ai.  (2001),  show  that  the  time 
response  of  PMS  imaging  probes  is  limited  by  the  time  constant  of  the 
photodiode  array  and  first  stage  amplifier,  leading  to  a  significant  reduction  in 
sensitivity  to  particles  with  sizes  <  -  75  pm  at  airspeeds  of  120  m  s'1  and  greater. 
Similar  results  have  been  obtained  using  a  high-speed  spinning  wire  to  test  the 
time  response  of  the  DMT  CIP  probe  (Nagel  2003  -  personal  communication)1 


2.1  2D-S  Probe  Description 

The  new  2D-S  probe  is  designed  to  minimize  errors  associated  with  DOF, 
sample  volume  uncertainties  and  time  response.  A  photograph  of  the  2D-S 
probe  illustrating  its  functional  aspects  is  shown  in  Figure  1.  Two  diode  laser 
beams  cross  at  right  angles  and  illuminate  two  linear  128-photodiode  arrays. 
The  laser  beams  precisely  define  a  rectangular  overlap  region  in  the  middle  of 
the  probe.  The  2D-S  simultaneously  records  two  independent  images  of  the 
same  particle  in  the  region  where  the  beams  overlap,  and  single  images  of  all 
particles  outside  the  overlap  region. 


1  Dr.  Dagmar  Nagel  of  the  Institut  for  Atmospharenphysik  (GKSS)  tested  CIP  probes  at  the 
Meteorology  Research  Institute  (MR1)  of  Japan  and  at  the  National  Center  for  Atmosphere 
Research  (NCAR)  and  found  that  at  75  m  s'1,  the  smallest  diameter  wire  that  could  be  imaged 
was  75  Mm. 
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Figure  1.  Photographs  of  the  SPEC  2D-S  probe.  See  text  for  description 


SPEC  2D-S  Probe 


DSP  Control  Board  128-Pfwtodiode  Naa tod  Sapphire 

with  Hardware  Data  Array  Board  with  Windows  with  KnUe- 
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and  Comparators 
(1  of  2  Boards) 


The  electro-optics  used  in  the  2D-S  are  on  the  cutting  edge  of  technology.  In 
order  to  substantially  increase  the  time  response  of  the  photodiode  array  and  first 
stage  amplifier,  a  new  custom  128-element  linear  photodiode  array  was  built  by 
CentroVision  Corporation.  Each  photodiode  in  the  array  is  50  pm  instead  of  200 
pm  used  in  the  conventional  (PMS  and  DMT)  photodiode  arrays,  thereby 
reducing  capacitance  and  increasing  speed.  The  128-photodiode  array  chip  is 
positioned  on  a  12-layer  printed  circuit  board  that  is  specially  designed  to  reduce 
lead  length  and  increase  speed  (Figure  1).  The  2D-S  probe  response  was 
tested  in  the  laboratory  and  compared  with  the  response  of  a  PMS  2D-C  probe. 
Figure  2  shows  digital  oscilloscope  traces  of  the  time  response  of  the  2D-C  and 
2D-S  probes  to  a  sharp  (5  ns  fall  time)  negative-going  pulse  of  laser  light  (i.e  , 
when  the  laser  beam  is  blocked  corresponding  to  when  a  particle  occults  a 
photodiode).  The  2D-C  has  a  (1/e)  response  time  constant  of  760  ns,  while  the 
2D~S  probe  has  a  response  of  41  ns,  which  is  nearly  20  times  faster  than  the  2D- 
C  probe. 

An  additional  test  that  demonstrates  the  extremely  fast  response  time  of  the 
2D-S  was  conducted  by  Dr  Dagnar  Nagel  with  a  unique  high-speed  spinning 
wire  calibration  unit  developed  and  operated  by  the  Institut  fur 
Atmospharenphysik  (GKSS).  Using  the  GKSS  calibration  unit,  the  2D-S  was 
able  to  completely  image  an  8  wire  spinning  at  a  velocity  of  233  m  s'1. 
Figure  3  shows  a  photograph  of  the  GKSS  spinning  wire  set  up  and  Figure  4 
shows  a  digital  oscilloscope  recording  of  the  response  of  the  2D-S  to  the  8  pm 
wire  spinning  at  233  m  s'  .  In  comparison,  using  the  same  GKSS  spinning  wire 
set-up,  the  best  response  that  the  PMS  2D-C  and  DMT  CIP  probes  could 
achieve  was  to  image  a  75  )im  wire  at  75  m  s*1. 
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Figure  2.  Digital  oscilloscope  traces  and  computation  of  time  constants  from 
the  front  ends  of  (left)  PMS  2D-C  probe  and  (right)  SPEC  2D-S  probe.  The 
response  is  generated  by  illuminating  each  photodiode  array  with  pulsed  laser 
light;  the  falling  edge  simulates  shadowing  by  a  particle.  The  1/e  time  constant 
of  the  2D-S  is  41  ns,  nearly  20  times  faster  than  the  760  ns  time  constant  of 
the  2D-C  (note  that  the  time  scale  in  the  2D-C  plot  is  500  ns  per  box  and  it  is 
25  ns  per  box  on  the  2D-S  plot). 


Figure  3.  GKSS  high-speed  spinning  wire  tests  of  the  2D-S  probe. 
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Figure  4.  Oscilloscope  captured  electronics  waveforms  for  an  8  pm  wire 
spinning  through  the  sample  volume  at  233  m  s'1. 


A  12-layer  printed  circuit  board  with  a  32-bit  AD-2192  DSP  running  at  80  MHz 
controls  very  large,  ball-grid  array  Complex  Programmable  Logic  Devices 
(CPLDs)  that  perform  hardware  data  compression  to  handle  the  extremely  high 
data  rates  (Figure  1).  The  2D-S  probe  outputs  data  at  a  maximum  rate  of  20 
million  samples  per  second  per  128-photodiode  array.  This  amounts  to  a 
maximum  “raw"  data  rate  of  1 .28  Gbytes  s"1  Since  much  of  the  area  is  '“white 
space"  between  particle  images,  hardware  data  compression  reduces  this  rate  by 
about  a  factor  of  ten  to  a  thousand  (depending  on  particle  size  and 
concentration),  making  it  possible  to  transfer  data  via  a  high-speed  serial 
interface  (HDLC).  As  discussed  above,  the  standard  2D-C  imaging  probes  miss  a 
substantial  fraction  of  particles  <  75  pm  in  diameter.  In  the  2D-S,  the  a  nival  time 
of  each  particle  is  recorded  and  the  x-y  position  of  each  particle  in  the  overlap 
region  of  the  two  laser  beams  is  recorded.  Finally,  as  shown  in  Figure  1 ,  the  2D- 
S  probe  is  specially  designed  with  non-wetting,  heated  sapphire  windows  that 
have  a  knife-sharp  leading  edge  to  minimize  particle  bounce  and  breakup  that  is 
suspected  to  contaminate  imaging  probes  (Jensen  and  Granek  2002).  Data  from 
the  2D-S  are  recorded  through  standard  aircraft  PMS  2D-C  wiring  using  an 
Ethernet  connection  that  supports  data  transfer  rates  up  to  100  mbs 
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2.2  Results  of  Preliminary  Flight  Tests 


The  2D-S  probe  was  installed  on  the  SPEC  Learjet  and  flown  on  seven 
research  flights  during  November-December  2003  and  eight  flights  during  April- 
May  2004.  In  November  2004,  it  was  installed  on  the  National  Center  for 
Atmospheric  Research  (NCAR)  C-130  aircraft  for  the  Rain  In  Cumulus  over  the 
Ocean  (RICO)  project.  Data  were  collected  by  the  Learjet  in  cirrus  clouds,  wave 
clouds  and  upslope  clouds  along  the  Front  Range  of  the  Rocky  Mountains  in 
Colorado,  Data  were  collected  by  the  C-1 30  in  an  upslope  cloud  along  the  Front 
Range,  and  on  1 3  missions  in  warm  cumulus  clouds  near  Antigua  during  the 
RICO  project.  Figure  5  shows  photographs  of  the  2D-S  probe  installed  on  the 
Learjet  and  the  NCAR  C-130. 


Figure  5.  Photographs  of  the  2D-S  probe  installed  on  (left)  the  SPEC  Learjet  and 
(right)  NCAR  C-130. 


Here  we  present  comparisons  of  2D-S  images  with  images  collected  using  a 
cloud  particle  imager  and  conventional  imaging  probes  installed  on  the  research 
aircraft.  The  2D-S  results  are  based  on  preliminary  software  programs.  Because 
the  analysis  software  is  preliminary,  we  opt  to  show  only  particle  counts  in  all 
cases  except  one,  which  is  when  the  probe  was  in  small  drizzle  drops  during 
RICO.  In  the  RICO  case,  we  incorporated  recent  work  by  Korolev  (2005)  to 
correct  for  out-of-focus  images,  and  also,  because  the  drops  are  small  (<  1 50 
pm)  we  believe  that  artifacts  that  may  be  generated  from  drop  shattering  and/or 
incomplete  diffraction  rings  are  minimal.  The  development  of  mature  software 
capable  of  automatically  rejecting  artifacts  (e  g.,  images  fragmented  from  large 
particles  shattering  on  the  probe  tips),  accurate  determination  of  stereo  particles, 
sizing  corrections  to  non-spherical  particles  and  improving  image-to-volume 
algorithms,  is  a  work  in  progress. 
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2.2.1  Mixed-phase  Ups  lope  and  Wave  Clouds 

Figure  6  shows  examples  of  2D-S,  CPI  and  2D-C  particle  images  collected  at 
-13  0  C  by  the  NCAR  C-130  in  a  mixed-phase  upslope  cloud  east  of  the 
Colorado  Front  Range,  The  10  pm  resolution,  0.128-cm  width  of  the  sample 
area  and  high-speed  response  of  the  2D-S  probe  allow  it  to  show  ice  particles 
within  a  field  of  small  cloud  droplets.  Compared  to  the  2D-C  particle  images,  the 
2D-S  shows  considerably  more  detail  of  ice  crystal  structure.  Figure  6  also 
shows  examples  of  “stereo"  views  of  images  that  passed  through  the  overlap 
region  of  the  2D~S  probe.  The  stereo  views  reveal  details  of  the  three- 
dimensional  structure  of  particles  that  cannot  be  determined  from  a  single 
shadow-image  of  the  particle 
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Figure  6.  Examples  of  images  from  the  2D-S,  CPI  and  2D-C  probes  in  a  mixed- 
phase  upslope  cloud  investigated  by  the  NCAR  C-130  research  aircraft.  {All  of 
the  particle  images  in  this  figure  have  been  hand-picked  and  placed  using  a 
graphics  tool,  i.e.,  the  spatial  relationships  of  the  particles  are  not  as  observed  by 
the  probes.) 

Figure  7  shows  an  example  of  measurements  from  the  SPEC  Learjet  during 
penetration  of  a  wave  cloud  that  formed  over  the  Colorado  Front  Range.  The 
evolution  of  particles  is  shown  starting  from  the  top  of  the  figure,  which  is  the 
leading  (western)  edge  of  the  wave  cloud.  The  wave  cloud  undergoes  the 
expected  particle  evolution  at  -27  °C  (Baker  and  Lawson  2006),  i.e.,  a  transition 
from  supercooled  cloud  droplets  near  the  leading  edge,  changing  to  a  mixed- 
phase  region,  which  rapidly  turns  into  glaciated  cloud.  The  cloud  drop  size 
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distribution  from  an  FSSP  is  displayed  at  the  left  of  the  figure,  showing  that  the 
size  of  the  supercooled  drops  ranges  from  about  4  to  18  pm  in  diameter,  with  a 
mean  diameter  of  11.6  pm.  Both  the  CPI  and  the  2D-S  probes  detect  the  cloud 
drops  as  seen  in  the  sample  images.  The  location  of  the  first  ice  particles 
encountered  in  the  mixed-phase  cloud  region  is  of  interest  to  cloud  physicists, 
because  this  relates  to  ice  nucleation,  a  process  that  is  still  poorly  understood  in 
wave  clouds  (Baker  and  Lawson  2006).  The  data  in  Fig.  7  shows  that  the  2D-S 
detects  ice  particles  much  closer  to  the  cloud  leading  edge  than  either  the  CPI  or 
the  2D-C  probes.  The  2D-C  does  not  detect  particles  with  sizes  <  -  100  pm  due 
to  the  slow  response  of  the  photodiode  array  and  front-end  amplifier.  (Later  in 
this  section,  we  show  additional  probe  comparisons  from  the  RICO  project  that 
support  this  statement.)  The  CPI  easily  sees  the  cloud  drops,  but  since  the 
imaging  laser  is  triggered,  the  probe  response  is  dominated  by  the  high 
concentration  of  drops,  so  it  is  often  inactive  when  larger  (ice)  particles  enter  the 
sample  volume.  The  ability  of  the  2D-S  to  detect  small  ice  ‘'earlier"  into  the 
traverse  of  the  wave  cloud  may  provide  the  basis  for  a  better  understanding  of 
the  ice  nucleation  process  in  clouds. 


2.2.2  Cold  Wave  and  Cirrus  (Ice)  Clouds 

Figure  8  shows  examples  of  small  spheroidal  particles  and  budding  rosettes 
in  a  cold  (-  45  °  C)  cirrus  doud.  The  figure  is  configured  to  show  the  transition 
from  a  region  of  small  particles  to  budding  rosettes  and  highlights  how  the  2D-S 
probe  sees  these  particles  and  the  2D-C  barely  responds.  This  is  emphasized 
by  the  total  partide  count,  where  the  2D-S  recorded  268,897  partides  compared 
with  264  partides  seen  by  the  2D-C  over  a  40  s  time  period.  Also,  the  shapes  of 
many  of  the  larger  (>  -  100  pm)  budding  rosettes,  which  are  confirmed  by  the 
CPI  images,  are  discemable  and  can  be  identified  as  rosettes  in  the  2D-S 
imagery,  whereas  the  2D-C  sees  only  1-5  pixel  images  whose  shape  is 
indistinguishable.  Figure  9  shows  another  example  of  small  ice  partides  and 
budding  rosettes  in  a  cold  (-  50  0  C)  drrus  doud.  In  this  figure,  six  2D-S  image 
strips  have  been  aligned  horizontally.  Notable  in  this  figure  is  the  large  number 
of  small  ice  partides  in  the  absence  of  large  partides  that  could  have  shattered 
on  the  probe  tips.  Observations  of  high  concentrations  of  small  ice  partides  in 
cirrus  have  been  reported  by  some  investigators  (Lawson  et  al.  2001 ;  Gayet  et 
al.  2002;  Lawson  et  al.  2006),  but  questioned  by  others  as  being  a  possible  result 
of  crystal  shattering.  The  2D-S  probe  may  shed  new  light  on  this  issue. 
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Figure  7.  Examples  of  2D-S,  CPI  and  2D-C  particle  images  observed  starting 
from  (top  of  figure)  the  leading  edge  of  a  wave  cloud  containing  supercooled 
drops,  extending  into  the  mixed-phase  region  and  terminating  (bottom  of  figure) 
in  the  glaciated  region.  The  FSSP  cloud  drop  distribution  near  the  leading  edge 
is  shown  at  the  left.  Examples  of  2D-S  stereo  images  are  shown  in  red. 
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Comparison  of  Images  of  Small  Ice  Particles 
from  2D-S,  CPI  and  2D-C  Probes 
in  a  Cirrus  Cloud 
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Figure  8.  Examples  of  images  from  the  2D-S,  CPI  and  2D-C  probes  in  a  anus 
cloud  (-45 0  C)  investigated  by  the  SPEC  Learjet. 
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Figure  9.  Example  of  images  of  budding  rosettes,  irregular-shaped  particles  and 
small  particles  recorded  by  the  2D-S  during  a  test  flight  in  a  thin  cirrus  cloud  at  - 
50  °C. 

Figure  10  shows  image  data  from  several  particle  imaging  probes  on  the 
SPEC  Learjet  and  WB-57F  research  aircraft,  operated  by  the  National  Aviation 
and  Space  Agency  (NASA).  Both  aircraft  participated  in  the  Mid-Latitude  Cirrus 
Experiment  (MidCiX)  and  flew  a  coordinated  mission  in  a  system  of  persistent, 
vertically  stacked  wave  clouds.  The  wave  cloud  system  was  located  over  the 
Colorado  Front  Range  and  had  a  cloud  base  temperature  of  -43 0  C  and  a  doud 
top  temperature  of -63  0  C.  A  Rosemount  icing  probe  (Mazin  et  al.  2001;  Cober 
et  al.  2001)  indicated  there  was  no  detectable  supercooled  liquid  water  drops  in 
these  douds. 
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Figure  10,  Comparison  of  particle  images  from  various  probes  installed  on  the 
SPEC  Learjet  and  NASA  WB-57  in  a  cold  wave  cloud  where  the  two  aircraft 
penetrated  the  same  region  of  cloud. 
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Figure  10  shows  portions  of  the  Lear  and  WB-57  flight  tracks  as  they 
penetrated  the  same  wave  doud,  based  on  video  records  and  flight  notes  from 
both  aircraft.  The  Learjet  penetrated  the  cloud  at  10.6  Km  (-54  °  C)  and  the  WB- 
57  penetrated  200  m  lower  in  the  doud  at  -53  0  C.  Both  aircraft  carried  CPIs 
and  2D-C  type  imaging  probes.  The  2D-C  probe  on  the  Learjet  was 
manufactured  by  PMS.  The  WB-57  flew  a  DMT  doud,  aerosol  and  predpitation 
spectrometer  (CAPS)  probe  (Baumgardner  et  al.  2001),  which  is  a  composite  of 
probes,  one  of  which  is  a  DMT  CIP.  The  CPI  partide  images  in  Fig.  10  show 
that  both  aircraft  penetrated  regions  with  small  spheroidal  partides  and  budding 
rosettes  (i.e.,  rosettes  with  short  branches).  The  maximum  dimension  of  the 
budding  rosettes  observed  by  the  Lear  is  about  75  pm,  while  the  WB-57  recorded 
some  slightly  larger  partides,  with  maximum  dimension  of  about  125  pm,  during 
its  penetration  200  m  lower  in  the  doud.  Of  note  is  the  observation  that  during 
the  entire  penetration  (and  adually  during  the  entire  flight),  there  were  no  partide 
images  recorded  by  the  Learjet  PMS  2D-C  probe.  The  DMT  CIP  probe,  on  the 
other  hand,  only  recorded  images  that  were  25  and  50  pm  (one  and  two  pixels). 
The  true  airspeed  of  each  aircraft  during  the  doud  penetrations  was  about  160 
ms'1.  Since  the  CPI  imaged  partides  up  to  125  pm  at  the  same  time  that  the 
CIP  recorded  25  and  50  pm  images,  it  would  appear  that  the  front-end  electro¬ 
optics  of  the  CIP  were  not  fast  enough  to  capture  the  entire  particle  image. 
Figure  10  shows  that  the  2D-S  probe  recorded  partide  images  as  small  as  10 
pm  and  with  tie  largest  image  dimension  being  200  pm  during  the  penetration. 
However,  since  the  2D-S  electro-optics  are  fast  enough  to  accurately  capture  the 
image,  it  is  possible  that  the  200  pm  images  are  out-of-focus  partides  that  are 
actually  about  125  pm  in  size  (see  Korolev  et  al.  1998). 

2.2.3  Warm  Cumulus  (Water)  Clouds 

The  RICO  project  addressed  several  aspeds  of  the  warm  rain  process,  one 
of  which  is  evolution  of  the  droplet  spectra  in  small,  warm  cumulus  douds 
(Rauber  et  al  2003).  The  2D-S  probe  can  detect  and  size  drops  from  10  pm  to 
>  1  mm  in  diameter,  so  it  is  capable  of  measuring  the  evolution  of  the  drop  size 
distribution  in  the  critical  size  region  where  doud  drops  coalesce  into  drizzle  and 
eventually  become  raindrops.  Measurements  of  the  drop  size  distribution  in 
warm  cumulus  have  been  made  during  several  previous  projects  using  PMS  2D- 
C  and  260X  probes  (e.g.,  Blyth  et  al  2003;  Jensen  and  Granek  2002).  We  show 
here  that  these  instruments  appear  to  have  improperly  measured  the  critical  size 
range  of  drops  with  diameters  from  about  25  to  150  pm. 

Figure  11  shows  examples  of  2D-S  and  2D-C  images,  along  with  PSDs  from 
the  FSSP,  2D-S  and  2D-C  during  a  250-s  penetration  in  a  doud  region  on  7 
January  2005.  During  this  time  period,  the  C-130  encountered  doud  and  drizzle 
drops  that  produced  in-focus  2D-S  images  that  extended  out  to  150  pm  in 
diameter.  The  FSSP  and  2D-C  data  were  processed  using  standard  NCAR 
processing  software  (NCAR  2003).  The  2D-C  depth-of-field  correction  for 
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sample  volume  comes  from  the  PMS  2D-C  manual  (NCAR  2003),  which  uses  a 
value  of  c  »  6  in  (1).  The  2D-S  images  are  shown  exactly  as  they  were  recorded. 
Based  on  theory  shown  in  Korolev  et  al.  (1998)  and  algorithmic  simulations 
recently  developed  for  SPEC  by  Korolev  (2005),  out-of-focus  2D-S  images  (i.e., 
“donuts’)  with  maximum  dimensions  £  30  pm  were  automatically  resized.  The 
Korolev  (2005)  algorithm  uses  the  ratio  Darea  /  Dma*,  where  Dafea  and  D™*  are 
explained  in  detail  in  Korolev  et  al.  <1998). 2  The  Korolev  (2005)  algorithm  is  only 
applied  to  the  portion  of  the  2D-S  PSD  &  30  pm  because  smaller  images  do  not 
contain  enough  information  (e.g.,  Poisson  spot)  to  reconstruct  the  original  particle 
size.  The  concentration  of  2D-S  images  <  20  pm  was  adjusted  using  depth-of- 
field  and  sample  volume  corrections  based  on  (1)  with  a  value  of  c  =  8,  which 
was  determined  from  laboratory  calibrations  with  glass  beads.  No  attempt  was 
made  to  eliminate  artifacts,  however,  in  the  RICO  PSD  presented  here,  there  is 
little  evidence  of  drop  shattering  or  other  artifacts.  Also,  the  occurrence  of 
images  “joined  together  is  minimal  and  can  be  ignored. 


Total  2D-S  Images  Recorded  =  314758 


Figure  11.  Comparisons  of  (left)  PSDs  from  FSSP,  2D-S  and  2D-C  and  (right) 
droplet  images  for  the  time  period  from  185340  to  185750  on  7  January  2005. 


2D-C  images  from  50  to  1 50  pm  (2  to  6  pixels)  cannot  be  resized  due  to 
the  poorer  instalment  response  time  and  pixel  resolution.  The  2D-C  drop  size 
distribution  was  computed  using  depth-of-field  corrections  published  in  the  PMS 
2D-C  manual.  As  shown  in  the  Table  1,  the  2D-C  recorded  a  total  of  46 
images  while  the  2D-S  horizontal  channel  recorded  314,758  images  during  the 
same  250-s  time  period;  90  of  the  2D-S  images  had  diameters  of  100  pm  or 
greater.  A  comparison  of  the  number  of  particles  in  Fig.  11  shows  that  the  2D-C 
recorded  about  0.1%  of  the  particles  seen  by  the  2D-S  in  the  size  range  from 


2  The  algorithm  developed  by  Korolev  is  described  in  Korolev  (2005),  which  can  be  viewed  at 
www.specmc.com.  and  will  be  the  subject  of  a  future  paper  by  Korolev  and  co-authors. 
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about  30  to  150  jim,  and  the  2D-C  PSD  was  about  an  order  of  magnitude  lower 
than  the  2D-S  in  this  size  range. 


Table  1.  Comparison  showing  response  of  2D-S  probe  horizontal  (H-channel) 
with  expected  counts  from  2D-C  probe  assuming  theoretical  instantaneous  time 
response  and  actual  counts  for  time  period  with  maximum  drop  diameter  of  (left) 
150  pm  and  (right)  maximum  drop  diameter  of  120  pm.  Data  shown  at  left  is 
from  RICO  for  250-s  time  period  from  185340  to  185750  on  7  January  2004,  and 
data  shown  at  right  is  a  continuous  65-s  time  period  within  the  250-s  time  period. 
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Of  the  total  images  recorded  by  the  2D-S,  30  of  the  images  had  diameters  > 
100  pm  and  the  2D-C  probe  did  not  record  any  images  >  100  pm.  Therefore, 
even  though  a  single  channel  of  the  2D-S  has  1.75  times  the  sample  volume  of 
the  2D-C  for  drops  <  150  pm,  the  absence  of  recorded  2D-C  images  >  100  pm 
cannot  be  explained  by  sampling  statistics  (Table  1).  The  260X  probe  recorded 
only  6  drops  during  the  same  250-s  time  period. 


The  2D-C  probe  did  not  record  any  images  larger  than  100  pm  in  diameter 
during  the  250-s  period,  which  suggests  that  the  few  images  it  did  record  may 
actually  have  been  out-of-focus  images  of  larger  (100  to  150  pm)  drops,  This  is 
supported  by  the  data  shown  in  Table  1,  which  shows  the  2D-S  horizontal  (H-) 
channel  counts  per  bin  for  65-s  during  the  250-s  period  when  the  2D-C  did  not 
record  any  images.  During  this  65-s  period,  the  2D-S  probe  H-channel  recorded 
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55,544  drop  images,  375  with  diameters  from  60  to  120  pm  and  12  with 
diameters  100  pm  and  larger.  Since  the  2D-C  probe  did  not  record  any  images, 
this  suggests  that  the  threshold  drop  size  that  the  2D-C  probe  detected  was 
about  >  100  pm  at  the  103  m  s'1  airspeed  of  the  C-130.  With  this  information, 
we  can  hypothesize  that  the  images  recorded  by  the  2D-C  in  the  50  to  100  pm 
size  bins  during  the  250-s  time  period  (Fig.  11)  are  actually  drops  in  the  size 
range  from  100  to  150  pm  that  produce  smaller  images  due  to  inadequate  time 
response  of  the  probe. 

Another  interesting  observation  drawn  from  Fig.  11  is  that  if  one  extrapolates 
the  PSD  from  the  FSSP  to  the  2D-C,  they  form  a  relatively  straight  line,  implying 
that  the  2D-C  PSD  is  reasonable.  Given  that  it  is  likely  that  the  2D-C  is  recording 
measurements  in  the  25  to  100  pm  size  range  that  are  actually  out-of-focus 
images  of  larger  drops,  one  can  draw  the  conclusion  that  many  of  the  PSDs 
previously  reported  in  this  size  range  are  erroneous.  To  be  more  specific, 
measurements  from  the  2D-C  with  time  constants  similar  to  the  NCAR 
instrument,  when  operated  at  100  m  s'1  or  faster,  likely  do  not  detect  (spherical 
water)  particles  with  sizes  <  100  pm.  If  particles  >  100  pm  are  present,  the  out- 
of-focus  images  from  these  particles  may  artificially  create  a  pseudo  PSD  from 
25  to  100  pm. 

A  comparison  of  measurements  from  the  (independent)  vertical  and 
horizontal  2D-S  channels  is  similar  to  comparing  the  response  to  two 
independent  imaging  probes.  Previous  efforts  to  compare  the  output  from 
several  2D-C  imaging  probes  installed  on  the  same  aircraft  have  shown  that 
errors  in  concentration  of  up  to  50%  are  expected  (Gayet  et  al.  1993).  Drop  size 
distributions  and  time-series  of  drop  concentration  from  both  the  vertical  and 
horizontal  2D-S  imaging  channels  are  shown  in  Fig.  12,  The  good  agreement  in 
the  time-series  of  drop  concentrations  between  the  vertical  and  horizontal 
channels  is  encouraging.  This  level  of  agreement  is  not  always  present  in  the 
RICO  data  set,  because  the  functioning  of  one  channel  or  the  other  is  sometimes 
degraded  in  the  salty  tropical  environment  due  to  window  contaminants,  and/or 
optical  alignment  issues. 
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Figure  12.  Comparison  of  (top)  time-series  data  and  (bottom)  scatter  plot  for  the 
“vertical"  and  “horizontal”  channels  of  the  2D-S  probe.  Occasional  data  dropouts 
in  the  time  series  are  when  a  channel  is  in  overload.  Each  square  in  the  scatter 
plot  is  a  4-s  average.  Measurements  are  for  the  cloud  penetration  from  185340  to 
185750  on  7  January  2005. 
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3.  Use  of  the  2D-S  on  NASA  Field  Campaigns. 

The  Navy  gave  permission  for  SPEC  to  install  the  2D-S  on  the  NASA  DC-8 
for  the  NASA  African  Monsoon  Multidisciplinary  Analysis  (NAMMA)  project,  and 
on  the  NASA  WB-57F  for  the  Costa  Rica  -  Aura  Validation  Experiment  (CR¬ 
AVE).  The  2D-S  performed  exceptionally  well  in  both  of  these  projects.  A  peer- 
reviewed  paper  that  discusses  data  collected  by  the  2D-S  installed  on  the  NASA 
WB-57F  in  subvisible  cirrus  clouds  in  the  tropical  tropopause  layer  is  contained  in 
Appendix  A. 

4.  Development  of  the  3V-CPI 

As  explained  in  Section  1,  the  Navy  3V-CPI  is  being  developed  in  parallel 
with  the  NSF/NCAR  HIAPER  3V-CPI.  A  Critical  Design  Review  (CDR)  was 
presented  to  the  NCAR  in  June  2006.  The  CDR  is  equally  applicable  to  the 
development  of  the  Navy  CPI  and  is  reproduced  below.  Since  the  time  of  the 
CDR  until  the  termination  of  the  Phase  II  contract  (2  July  2007),  development  of 
the  3V-CPI  has  continued  and  all  mechanical,  electrical  and  optical  components 
have  been  designed  and  a  large  majority  of  the  parts  were  fabricated.  Since  2 
July  2007,  SPEC  has  continued  development  of  the  Navy  3V-CPI  with  internal 
research  and  development  funds.  Initial  laboratory  tests  of  the  subsystems  are 
proceeding  and  final  tests,  including  flight  tests  on  the  HIAPER  G-V  are  expected 
to  be  completed  by  June  2008,  at  which  time  SPEC  plans  delivery  of  the  3V-CPI 
to  the  Navy  for  application  in  the  VOCALS-REx  project  planned  for  October 
2008. 


3V-CPI  Critical  Design  Review 
(Presented  to  the  Research  Aviation  Facility  at  the 
National  Center  for  Atmospheric  Research  in  June  2006): 
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Solidworks  Drawing  Showing  Internal 
Structure  of  2D-S/CPI  with  2D-S  Electro- 
Optics  and  Basler  504K  Digital  Camera 
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functions,  we  expect  a  slight  reduction  in  maximum 
image  throughput  per  second;  however  the  achieved 
throughput  exceeds  the  camera  capabilities 
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FIBER  /  POWER  INTERFACE  BOX  DETAILS 
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CONFLICT  BETWEEN  FIBER/POWER  INTERFACE  BOX  AND  3V-CPI 
CONNECTORS  WITHIN  NCAR  PROPOSED  MOUNT  STRUCTURE 


MOUNT  STRUCTURE  NEEDS  TO  ACCOMMODATE  INSTRUMENT  CONNECTORS  AND 
FIBER  /  POWER  INTERFACE  BOX-  CABLE  PASS  TO  DOGHOUSE  /  CABIN 
NCAR  PRELIMINARY  MOUNT  STRUCTURE  NEEDS  TO  BE  MODIFIED 
ACCESS  REQUIRED  TO  MATE  CONNECTORS  AND  REACH  MOUNTING  HOLES  IN 
BOTTOM  OF  PYLON 
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DIFFICULT  TO  INTERPRET  IMAGES  /  OUT  OF  FOCUS  EFFECTS, 
2D-S  SOFTWARE  IS  NOT  DIRECTLY  APPLICABLE 


NEW  OPTICAL  SYSTEM  ORIENTATION 
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NEW  ORIENTATION  EFFECT  ON  CPI  IMAGE  FOCUS 
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2D-S  OPTICAL  SYSTEM  PARAMETERS 
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TRIGGERING  ELECTRONICS  DESIGNED  TO  COMPENSATE  FOR  SOME 
THERMAL  SHIFT  IN  45  AND  90  OBJECT  PLANES 


2D-S  IMAGING  OPTICAL  SYSTEM,  EXISTING  VS  PROPOSED 
VS  NEWLY  PROPOSED  4  LENS  SYSTEM 


TOTAL  TRACK  LENGTH  =276  MM 


NEWLY  PROPOSED  2D-S  4  LENS  SYSTEM 
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3V-CPI  OPTO-MECHANICAL/PACKAGING  CONCERNS 
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DUAL  FOLD  4  LENS  SYSTEM  PROTOTYPED  IN  LAB 
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3V-CPI  90°  2D-S  SYSTEM 
ASSEMBLY  LAYOUT 
(45°  ALMOST  IDENTICAL) 


3V-CPI  COMBINED  OPTICAL  PATHS 


CPI  IMAGE  SYSTEM  COMPARISON,  M=5X 
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CPI  IMAGE  SYSTEM  LAYOUT  INCORPORATING  135  0  FOLD  PATH 


PTICAL  SYSTEM 
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System  performance  evaluated  -  with  water  drops,  glass  beads,  moving  wire. 
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EXAMPLE  OF  2D-S  LASER  BEAM  SHAPING 
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Schematics  delivered  as  part  of  this  design  review 
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AIRCRAFT  MEASUREMENTS  OF  MICROPHYSICAL  PROPERTIES  OF 
SUBVISIBLE  CIRRUS  IN  THE  TROPICAL  TROPOPAUSE  LAYER 
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Abstract.  Subvisible  cirrus  (SVC)  clouds  are 
often  observed  within  the  tropical  tropopause  layer 
(TTL)  and  have  been  shown  to  have  a  significant 
impact  on  the  earth  radiation  budget.  The  Costa 
Rica  Aura  Validation  Experiment  (CR-AVE) 
sponsored  by  the  National  Aeronautics  and  Space 
Administration  (NASA)  took  place  near  San  Jose, 
Costa  Rica  from  14  January  -  15  February  2006. 
The  NASA  WB-57F  sampled  SVC  in  the  TTL  from 
-75°C  to  -90°C  with  an  improved  set  of  cloud 
particle  probes.  The  first  digital  images  of  ice 
particles  in  the  TTL  are  compared  with  replicator 
images  of  ice  particles  collected  in  1973  by  a  WB- 
57F  in  the  TTL  The  newer  measurements  reveal 
larger  particles,  on  the  order  of  100  pm  compared 
with  <  50  pm  from  the  earlier  measurements,  and 
also  different  particle  shapes.  The  1973  particles 
were  mainly  columnar  and  trigonal,  whereas  the 
newer  measurements  are  quasl-spherical  and 
hexagonal  plates.  The  WB-57F  also  measured 
very  high  water  vapor  contents  with  some 
instruments,  up  to  4  ppmv,  and  aerosols  with 
mixed  organics  and  sulfates.  It  is  unknown 
whether  these  ambient  conditions  were  present  in 
the  1973  studies,  and  whether  such  conditions 
have  an  influence  on  particle  shape  and  the 
development  of  the  large  particles.  A  companion 
paper  (Jensen  et  al.  2007}  presents  crystal  growth 
calculations  that  suggest  that  the  high  water  vapor 
measurements  are  required  to  grow  ice  particles 
to  the  observed  sizes  of  100  pm  and  larger. 


1  Introduction 

Subvisible  cirrus  (SVC)  clouds  appear  mostly  in 
the  tropical  tropopause  iayer  (TTL)  up  to  the 
tropopause  (Wang  et  al.  1996;  Beyerle  et  al.  1998; 
McFarquhar  et  al.  2000;  Comstock  et  al.  2002; 
Clark  2005),  but  have  also  been  observed  at 
midlatftudes  (Sassen  and  Campbell  2001;  Immler 
and  Schrems  2002).  The  TTL  is  the  region  of  the 
tropical  atmosphere  that  lies  between  the  top  of 
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the  main  cumulus  outflow  layer  and  the  thermal 
tropopause  (Folkins  et  al.  1999).  This  layer  is  a 
transition  layer  between  dynamical  control  of  the 
vertical  mass  flux  by  tropospheric  convection,  and 
by  the  stratospheric  Brewer-Dobson  circulation, 
and  is  crucial  to  understanding  the  dehydration  of 
air  entering  the  stratosphere. 

Although  the  radiative  forcing  of  optically  thin 
SVC  is  relatively  small,  the  clouds  generally  cover 
a  large  horizontal  extent  and  are  considered  to  be 
radiatively  significant  (McFarquhar  et  al.  2000; 
Comstock  et  al.  2002).  Using  a  nadir-pointing 
Nd;YAG  lidar  on  a  Learjet,  McFarquhar  et  al. 
(2000)  report  that  a  spatially  thin  layer  of  cirrus, 
with  both  base  and  top  above  15  km,  was 
observed  in  the  central  Pacific  Tropics  29%  of  the 
time,  with  a  mean  thickness  of  0.47  km.  Based 
on  vertical  profiles  and  a  radiative  transfer  model, 
they  calculate  observed  heating  rates  of  up  to  1.0 
K  day'1,  principally  in  the  infrared,  and  cloud 
radiative  forcing  of  up  to  1.2  W  m'2  SVC  clouds 
are  occasionally  observed  above  the  tropopause 
at  midlatitudes  (Goldfarb  et  al.  2001)  and  in  Polar 
Regions  (Lelieveld  et  al.,  1999;  Karcher  and 
Solomon,  1999). 

Supersaturation  with  respect  to  ice  has  been 
observed  in  the  upper  troposphere  and  lowermost 
stratosphere  (Murphy  et  al.  1990;  Gierens  et  al. 
1999;  Ovartez  et  al.  2000;  Jensen  et  al.  1999, 
2001,  2005;  Peter  et  al.  2006).  Brewer  (1949) 
showed  that  air  entering  the  stratosphere  is  freeze 
dried  as  It  crosses  the  cold  tropical  tropopause. 
Jensen  et  al.  (1996)  suggested  that  SVC  can 
effectively  freeze-dry  air  crossing  the  tropical 
tropopause  to  the  observed  lower  stratospheric 
water  vapor  concentrations,  which  may 
significantly  affect  the  Earth's  radiation  budget 
(Comstock  et  al,  2002). 

SVC  clouds  were  studied  with  the  NASA  WB- 
57F  research  aircraft  in  January  -  February  2006 
during  the  Costa  Rica  -  Aura  Validation 
Experiment  (CR-AVE 

htto://cloud1.  arc,  nasa.gov/ave-costarica21. 

Figure  1  shows  a  photograph  of  the  WB-57F  and 
an  example  of  an  SVC  iayer  as  seen  from  the 
cockpit  of  the  aircraft.  While  there  have  been 
several  studies  of  SVC  using  airborne,  satellite 
and  ground-based  Ildar,  there  has  been  a  relative 
dearth  of  in  situ  microphysical  measurements, 
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Heymsfield  (1986)  and  McFarquhar  et  al.  (2000) 
discuss  measurements  made  by  a  WB-57F  in  a 
SVC  cloud  at  -84°C  over  the  Marshall  Islands  on 
17  December  1973,  Based  on  an  examination  of 
replicator  data,  they  concluded  that  particle  sizes 
ranged  from  about  5  to  50  pm  and  particle  shapes 
were  mostly  columnar  and  trigonal  plates,  with 
some  hexagonal  plates.  The  only  other 
microphysical  measurements  in  SVC  that  we  are 
aware  of  have  been  made  with  scattering  probes, 
such  as  the  forward  scattering  spectrometer  probe 
(FSSP)  For  example,  Thomas  et  al.  (2002)  found 
a  modal  peak  in  the  size  distribution  between 
about  2  and  17  pm,  and  Peter  et  al.  (2003)  found  a 
model  peak  of  10  to  12  pm,  both  studies  used  an 
FSSP-300  (Baumgardner  et  al.  1992),  which 
measured  particles  from  0.4  to  23  urn  In  size. 


Fig  1.  Photograph  of  (top)  Subvisible  Cirrus  (SVC) 
layer  taken  from  the  cockpit  of  the  NASA  WB-57F 
(shown  at  bottom  before  sensors  were  installed) 
during  the  transit  flight  from  Houston,  Texas  for 
the  Costa  Rica  Aura  Validation  Experiment  (CR¬ 
AVE).  See  acknowledgements  for  photo  credit. 

The  microphysical  measurements  collected  in 
SVC  by  the  WB-57F  research  platform  include,  for 
the  first  time,  digital  images  of  particles  horn  a 
cloud  particle  imager  (CPI)  probe  (Lawson  et  al. 
2001}  a  cloud  and  aerosol  particle  spectrometer 
(CAPS)  probe  (Baumgardner  et  al.  2001),  and  a 
new  device,  the  2D-S  probe  (Lawson  et  al.  2006). 


Data  from  the  CPI,  CAPS  and  2D-S  are  analyzed 
from  four  WB-57F  missions  when  all  three 
instruments  were  functioning.  Temperature, 
pressure  and  true  airspeed  measurements  come 
from  the  WB-57F  meteorological  measurement 
system  (Scott  et  al.  1990)  and  water  vapor 
measurements  are  from  the  Harvard  Integrated 
Cavity  Output  Spectroscopy  (ICOS)  instrument 
(Sayres  2006).  The  data  presented  here 
represent  approximately  1,800  km  of  WB-57  F 
flight  in  SVC  during  the  CR-AVE  project.  The  CPI 
recorded  8524  images  and  the  2D-S  probe 
recorded  over  100,000  images  of  ice  particles  in 
SVC. 

2  Instrumentation  and  data 

A  highlight  of  CR-AVE  field  campaign  was  the  full 
complement  of  state-of-the  art  microphysica) 
probes  used  to  extensively  investigate  SVC,  The 
combination  of  particle  probes,  CPI,  2D-S  and 
CAPS,  provided  extensive  overlap  in  particle  size 
distributions  and  high-resolution  digital  images  of 
the  particles  in  SVC.  Data  for  this  paper  were 
extracted  from  the  CR-AVE  archives.  Before 
being  placed  in  the  archives,  2D-S  data  were 
averaged  over  1  s,  CPI  data  were  averaged  over 
10  s  and  CAPS  data  were  averaged  over  30  s.  In 
order  to  generate  1-Hz  CAPS  data  that  could  be 
combined  with  2D-S  data,  each  30-s  averaged 
CAPS  value  was  repeated  30  times.  Both  2D-S 
and  CAPS  data  were  then  used  to  generate  a 
combined  particle  size  distribution  (PSD)  and 
derived  parameters  on  a  1  Hz  basis.  Since  the 
ice  particles  did  not  exceed  165  pm  in  this  study, 
particle  shattering  on  probe  inlets  is  not 
considered  to  be  an  issue  (Field  et  al.  2003). 

The  CAPS  contains  two  particle  probes,  a 
cloud  and  aerosol  spectrometer  (CAS),  which 
measures  forward  and  backward  scattered  light  to 
count  and  size  particles  from  0.5  to  50  pm,  and  a 
cloud  imaging  probe  that  images  particles  from  75 
to  1600  pm  with  25  pm  pixels.  The  CPI  has  2.3 
pm  pixels,  which  provides  high-resolution  digital 
images  of  ice  partides  in  SVC,  but  the  relatively 
small  CPI  sample  volume  (a  maximum  of  0.7  L  s'1) 
compared  with  optical  array  probes,  limits  its 
ability  to  adequately  represent  the  entire  particle 
size  distribution.  The  2D-S  is  a  new  imaging 
probe  that  has  true  10  pm  per  pixel  resolution  at 
jet  aircraft  airspeeds,  whereas  the  response  of  the 
older  2D  imaging  probes  Is  degraded  with 
airspeed.  At  the  true  airspeed  of  the  WB-57F 
(about  180  m  s'1  at  55,000  ft),  the  older  2D 
imaging  probes  do  not  detect  particles  with  sizes 
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less  than  50  to  100  (Lawson  et  al.  2006). 
Figure  2  shows  photographs  of  the  CPI  and  2D-S 
probes  that  were  installed  on  the  WB-57F  for  CR¬ 
AVE. 


Fig  2.  Photographs  of  the  CPI  and  2D-S  probes 
installed  under  the  left  wing  of  the  NASA  WB-57F 
for  the  CR-AVE  project. 

The  true  10  pm  resolution  of  the 
2D-S  facilitates  resizing  of  out-of-focus  images, 
which  has  otherwise  been  a  source  of  error  that 
has  hampered  older  2D  imaging  probes  (Korolev 
et  al.  1998),  Examples  of  2D-S  images  in  SVC 
and  re-sizing  using  an  algorithm  developed  by 
Korolev  (2007)  are  shown  in  Fig.  3.  The  algorithm 
is  based  on  theoretical  diffraction  and  has  been 
evaluated  using  images  of  glass  beads  in  the 
laboratory.  The  Korolev  algorithm  was  tested  on 
several  bead  sizes,  and  it  performed  better  on 
some  sizes  than  it  did  on  others,  for  reasons  that 
are  unknown  at  this  time.  Later  in  this  paper  (Fig. 
10)  we  show  results  of  a  PSD  that  has  been 
resized  using  the  Korolev  (2007)  algorithm 
compared  with  a  PSD  computed  using  only  in¬ 
focus  2D-S  images. 


Water  vapor  was  measured  with  multiple 
instruments  on  the  WB-57F  during  the  CRAVE 
flights.  The  Harvard  water  vapor  instrument 
photodissociates  H20  molecules  with  121.6  nm 
(Lyman-a)  radiation  and  detects  the  resulting  OH 
photofragment  fluorescence  at  315  nm  (Weinstock 
et  al.,  1994).  The  Harvard  ICOS  instrument  uses 
a  mid-infrared  quantum  cascade  laser  at  1484  cm 
1  to  obtain  measurements  of  H20,  as  well  water 
vapor  isotopes,  with  an  effective  path  length  of 
nearly  4  km  (Sayres  2006).  The  Jet  Propulsion 
Laboratory  Laser  Hygrometer  (JLH)  is  a  tunable- 
diode  laser  focused  that  operates  by  harmonic 
wavelength  modulation  spectroscopy  (May  1998). 
For  over  50  years  measurements  of  water  vapor 
concentration  in  the  TTL  have  been  conducted 
using  balloon-borne  instrumentation,  typically 
frost-point  chilled  mirrors  (Oilmans  et.  al.  2000).  A 
long-standing  discrepancy  has  persisted  between 
the  frostpoint  H20  measurements  and  aircraft- 
borne  H20  measurements  made  using  a  variety  of 
techniques.  This  discrepancy  was  particularly 
glaring  under  the  very  dry  tropopause  conditions 
sampled  during  CR-AVE,  with  the  frostpoint  and 
JLH  measurements  typically  reporting  1-2  ppmv 
whereas  the  ICOS  and  Harvard  Water  Vapor 
instruments  report  3-4  ppmv.  This  discrepancy  in 
water  vapor  measurements  is  germane  to  this 
paper,  because  calculations  shown  in  Jensen  et 
al.  (2007)  suggest  that  water  vapor  limits  the 
maximum  size  to  which  ice  particles  can  grow 
before  they  fall  out  of  the  upper  region  of  the  TTL. 
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Fig  3.  2D-S  particle  images  from  the  vertical  (V) 
channel  of  the  2D-S  probe  and  an  example  of  size 
correction  of  a  donut  image  using  the  Korolev 
(2007)  resizing  scheme 
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3  Characteristics  of  subvisible  cirrus 

3.1  Spatial  and  vertical  extent 

Figure  4  shows  the  WB-57F  flight  track  on  the 
transit  flight  (TF  1)  to  San  Jose,  Costa  Rica  on  14 
January  2006,  and  flight  tracks  for  five  research 
flights:  RF  1  (17  January  2006),  RF  5  (25  January 
2006),  RF  8  (1  February  2006),  RF  9  (2  February 
2006)  and  RF  11  (7  February  2006).  These  are 
flights  when  there  were  CPI  data  available.  CPI 
data  were  used  to  identify  regions  when  the  WB- 
57F  was  in  SVC,  because  out  of  ail  the  particle 
probes,  it  was  the  most  reliable  indicator  of  SVC. 
CPI  data  were  also  used  to  determine  the  shape 
of  ice  particles,  because  of  its  high-resolution 
imagery.  Both  CAPS  and  2D-S  data  were 
available  for  TF  1,  RF  5,  RF  8  and  RF  9,  and  data 
from  these  missions  were  used  to  compute  the 
mean  PSD's  and  average  bulk  parameters 
computed  in  SVC.  Also  shown  in  Fig.  4  are 
vertical  profiles  of  the  transit  flight  and  research 
flights  where  the  CPI  observed  SVC.  The  data  in 
Fig.  4  show  that  SVC  was  frequently  observed  and 
that  it  was  consistently  found  between  about 
53,000  ft  (16.2  km)  and  58,000  ft  (17.7  km)  msl, 
which  corresponds  to  a  temperature  range  from 
about  -  75  to  -  85°C. 

As  suggested  in  Fig.  4,  the  WB-57F  was 
generally  in  and  out  of  SVC  clouds  during  its 
mission  profile,  and  sometimes  porpoised  from  the 
top  to  the  bottom  of  the  cloud  layer.  To  avoid 
averaging  regions  of  clear  air  with  SVC  cloud,  the 
missions  were  divided  into  flight  legs  that  were 
chosen  when  the  average  2D-S  particle 
concentration  was  >  5  L'1  for  5-km  or  greater 
without  a  containing  continuous  period  of  clear  air 
(2D-S  concentration  <  0  01  L'1)  that  was  1-km  or 
greater.  This  resulted  in  19  flight  legs  that  were 
about  25  to  200  km  in  length,  with  one  very  long 
leg  that  extended  for  2,400  km.  Figure  5  shows 
an  example  of  flight  legs  that  were  used  for 
making  PSD's  from  a  time  series  of  2D-S  particle 
concentration  during  a  portion  of  the  mission  on  2 
February  2006. 

3.2  Particle  habits 

Figure  6  shows  examples  of  CPI  images  in  SVC 
that  have  been  sorted  into  particles  <  65  pm  and 
particles  z  65  pm  in  maximum  dimension.  The  65 
pm  dimension  was  chosen  because  visual 
inspection  showed  that  the  shape  of  the  CPI 


images  typically  changed  at  about  this  size. 
Images  <  65  pm  appear  to  be  mostly  quasi- 
spherical  and  there  is  little  light  transmission 
through  the  images.  Images  a  65  pm  start  to 
appear  disk  shaped,  with  distinctive  edges  that  are 
not  apparent  in  the  smaller  spheroidal  images. 
Also,  several  of  the  images  £  65  pm  start  to  take 
on  a  plate-like  shape.  This  can  be  seen  more 
clearly  in  Fig.  7,  which  shows  expanded  views  of 
several  of  the  particles  and  reveals  the  distinctive 
edges  on  the  basal  faoes  of  the  particles,  dearly 
revealing  hexagonal  crystal  structure. 
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Fig  4.  (Top)  WB-57F  flight  tracks  during  CR-AVE 
showing  (in  blue)  regions  where  SVC  was 
encountered  and  (Bottom)  vertical  profile  showing 
location  of  SVC.  The  color  code  in  the  right  panel 
corresponds  with  the  research  flight  (RF)  number 
shown  in  the  legend. 
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Fig  5.  Example  of  a  time  series  of  WB-57F  2D-S 
total  particle  concentration  in  SVC  showing  4  of 
the  19  ‘legs'  {red  arrows)  used  in  the  analyses  in 
this  paper.  Legs  were  selected  using  an  objective 
criteria  explained  in  the  text. 

The  CPI  images  in  Figs.  6  and  7  are  in  sharp 
contrast  to  the  only  other  known  images  of  ice 
particles  in  the  TTL,  which  were  collected  on  17 
December  1973  using  a  replicator  installed  on  a 
WB-57F  research  aircraft  flying  between  16.2  and 
16.7  km  (-83°C  to  -84°C)  over  the  Marshall  Islands 
(Heymsfield  1986;  McFarquhar  et  al.  2000). 
Heymsfield  (1986)  shows  a  few  images  where  the 
shapes  of  the  ice  particles  are  mostly  columnar 
and  trigonal  (with  some  plates),  ranging  in  size 
from  about  5  to  50  pm. 

In  order  to  be  certain  that  the  CPI  images 
have  disk-like  edges  that  are  not  an  optical 
aberration,  we  captured  images  of  80  pm  glass 
beads  using  the  CR-AVE  CPI  in  the  laboratory  and 
examined  images  on  all  locations  in  the  CPI 
viewing  volume.  A  few  images  of  glass  beads  that 
come  the  closest  to  producing  what  appears  to  be 
an  edge  are  also  shown  in  Fig.  7.  It  is  apparent 
from  comparing  the  ice  particle  images  in  SVC 
with  the  images  of  spherical  glass  beads  that, 
most  likely,  the  edges  are  not  an  optical  aberration 
of  a  (glass)  sphere.  Also,  we  could  not  find  any 
evidence  of  disk-like  edges  upon  examining  quasi- 
spherical  images  of  50  to  100  pm  ice  particles  in 
anvil  cirrus  observed  during  CR-AVE,  further 
suggesting  that  the  disk-like  and  plate-like  SVC  ice 
particles  are  real. 

A  geometric  analysis  of  a  CPI  image 
(Appendix  A)  suggests  that  the  aspect  ratio  (basal 
to  prism  faces)  is  on  the  order  of  6:1.  However,  it 
is  not  possible  from  this  crystal  image  to  determine 
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if  the  aspect  ratio  is  actually  larger  than  6:1, 
because  of  the  amount  of  blur  in  the  image  due  to 
diffraction  and  optical  aberrations.  For  example, 
the  CPI  image  in  the  lower  left  frame  of  the  plate¬ 
like  images  in  Fig.  7  appears  to  have  an  aspect 
ratio  of  6:1,  assuming  it  is  a  plate  being  viewed 
orthogonal  to  the  prism  face  (and  1:6  if  it  is 
actually  a  column).  However,  if  it  is  a  tilted  plate,  it 
will  actually  have  a  slightly  larger  aspect  ratio 
Also,  if  the  prism  face  is  significantly  blurred  due  to 
diffraction  and  optical  aberrations,  the  dimension 
of  the  prism  face  will  be  artificially  lengthened,  and 
the  aspect  ratio  will  be  erroneously  decreased. 

The  aspect  ratio  of  these  larger  crystals  is 
germane  to  the  crystal  growth  calculations  and 
relative  humidity  measurements  presented  in  a 
companion  paper  (Jensen  et  al.  2007).  In  this 
paper  we  choose  not  to  offer  an  estimate  of  the 
largest  likely  aspect  ratio  of  these  crystals,  since 
any  estimate  will  contain  some  degree  of 
ambiguity  and  will  be  a  subjective  analysis  of 
thickness  of  the  image  edge. 

Figure  8  shows  histograms  of  particle  habits 
classified  by  a  trained  analyst  for  all  of  the  CPI 
images  collected  in  SVC  and  from  only  those 
images  >  65  pm.  The  histograms  comparison 
between  particles  <  65  pm  and  those  2  65  pm 
shows  considerably  fewer  quasi-spheroids  in  the 
larger  particle  histogram  and  a  significant  increase 
in  plate-like  and  edged  irregulars  (also  note  that 
the  “columns"  may  be  plates  viewed  on  edge). 

The  CR-AVE  data  in  SVC  observed  near 
Costa  Rica  contained  a  few  particles  that  were 
significantly  larger  (by  about  100  pm)  and 
markedly  different  crystal  habits  than  the 
observations  reported  by  Heymsfield  (1986)  over 
the  Marshall  Islands.  The  lack  of  trigonal  ice  in  the 
CR-AVE  SVC  data  set  is  curious.  The  replicator 
images  shown  in  Heymsfield  (1986)  are  clearly 
trigonal.  The  CPI  images  are  not  quite  as  sharp 
as  the  replicator  images,  but  once  particles  are 
larger  than  about  30  pm.  trigonal  shapes  will  be 
clearly  identifiable,  if  they  exist.  Thus,  the 
differences  in  particle  habits  between  the 
Heymsfield  (1986)  case  and  the  CR-AVE  data  set 
cannot  be  attributed  to  instrumental  effects.  The 
data  set  collected  in  1973  is  33  years  previous  to 
the  CR-AVE  data,  and  the  possibility  exists  that 
differences  in  the  moisture  field  and  atmospheric 
chemistry  have  changed  over  the  past  three 
decades;  however,  moisture  and  chemistry 
instrumentation  was  not  installed  on  the  WB-57F 
for  the  1973  flights. 
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Fig.  6.  Examples  of  CPI  images  in  subvisible  cirrus  (SVC)  of  (top)  ice  particles  >  65  pm  and  (bottom)  <.  65 


urn. 

3.3  Particle  size  distributions  (PSD’s)  and  bulk 
properties 

Figure  9  shows  PSD's  as  a  function  of  number, 
area  and  mass,  based  on  -1800  km  of  WB-57F 
flight  in  SVC  during  the  four  missions  discussed  in 
Sect.  3.1.  Compared  with  the  2D-S  mean  PSD, 
there  is  a  sharp  drop  in  the  CAPS  (CAS)  mean 
concentration  at  5  pm.  The  CAS  data  applies  Mie 


scattering  theory  to  size  particles,  and  implicit  in 
Mie  theory  is  the  assumption  that  the  particles  are 
perfect  spheres.  Since  the  CAS  is  sizing  ice  that 
is  not  perfectly  spherical,  the  sizing  uncertainty  is 
unknown.  Also,  it  uses  three  stages  of  gain 
amplifiers,  and  the  transition  from  one  of  the 
stages  to  another  occurs  at  5  pm  (where  the  dip 
starts  in  the  CAPS  PSD).  We  chose  to  use  the 
2D-S  PSD  in  the  size  region  from  10  to  30  pm 
region  and  interpolate  between  CAS  and  2D-S 
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between  5  and  10  pm.  Figure  9  shows  that, 
compared  with  the  2D-S,  the  CAPS  indicates  a 
higher  concentration  in  the  size  range  from  50  to 
150  pm,  and  that  it  recorded  particles  as  large  as 
225  pm,  while  the  largest  2D-S  particle  was  165 
pm.  Both  the  CAPS  and  2D-S  data  were 
processed  using  image-resizing  schemes  based 
on  Korolev  (2007).  A  check  on  2D-S  resizing  was 
accomplished  by  comparing  a  PSD  of  infocus 
images  with  resized  images  using  the  Korolev 
(2007)  technique.  However,  there  were  only  18 
2D-S  particles  >  100  pm  in  the  CR-AVE  data  set 
and  only  a  few  of  these  were  in  focus,  so  a 
comparison  using  the  CR-AVE  data  set  would  not 
produce  adequate  statistics.  Instead,  a  very  large 
data  set  (137,222  total  2D-S  images  with  5,601 
images  >  100  pm)  collected  in  drizzle  during  the 
2005  Rain  in  Cumulus  over  the  Ocean  (RICO) 
experiment  was  processed.  The  RICO  drizzle 
images  appear  very  similar  in  shape  to  the  quasL 
spherical  SVC  images. 

Figure  10  shows  the  results  of  comparing 
PSD's  using  only  infocus  2D-S  images  with  all 
images  processed  using  standard  2D  technique 
(i.e.,  the  maximum  image  dimension  across  the 
array)  and  images  resized  using  the  Korolev 
(2007)  scheme.  Appendix  B  discusses  sizing 
errors  using  only  infocus  images,  which  will 
usually  resuit  in  a  PSD  with  the  smallest  particle 
sizes.  In  Fig.  9,  using  infocus  images  does 
produce  the  smallest  PSD,  while  the  PSD 
produced  based  on  the  Korolev  (2007)  algorithm 
falls  between  the  Infocus  and  standard  processing 
PSD.  Based  on  the  comparisons  shown  in  Figs. 
9  and  10,  it  appears  that  the  2D-S  produces  more 
reliable  information  in  from  10  to  165  pm  (the 
largest  2D-S  particle  size  observed),  and  therefore 
the  2D-S  data  are  used  in  that  size  range. 

McFarquhar  et  al.  (2000)  show  data  in  SVC 
from  a  Particle  Measuring  Systems  (PMS)  axially 
scattering  spectrometer  probe  (ASSP)  with  a  size 
range  from  2  to  30  pm,  and  PMS  1-D  probe  that  is 
intended  to  size  particles  from  20  to  300  pm. 
However,  recent  studies  (Lawson  et  al.  2006) 
suggest  that  the  slow  time  response  of  the  1-D 
probe,  which  uses  a  diode  array  and  front-end 
amplifier  comparable  to  the  2D-C,  limit  its 
response  to  particles  larger  than  about  100  pm  at 
the  Learjet  airspeeds  used  in  that  study.  The  PSD 
shown  in  Fig.  1  in  McFarquhar  et  ai.  (2000) 
extends  from  2  to  15  pm,  with  some  1-D  probe 
activity  in  the  20  to  50  pm  size  range.  However, 
these  are  likely  to  be  much  larger  particles  that  are 
undersized  due  to  limitations  imposed  by  the  time 
response  of  the  1-D  probe.  Interestingly,  the 


ASSP  PSD  that  extends  from  2  to  15  pm  shown 
by  McFarquhar  et  al.  (2000)  is  similar  to  the 
composite  CR-AVE  PSD  shown  in  Fig.  9,  (but  it  is 
not  similar  to  the  CAS  PSD  alone).  Figure  9  also 
shows  average  composite  PSD's  as  a  function  of 
particle  area  and  particle  mass.  The  particle  area 
size  distribution  strongly  influences  optical  cloud 
properties,  such  as  extinction  and  optical  depth. 
The  particle  mass  size  distribution  has  an  impact 
on  ice  water  content  (IWC)  and  ice  water  path. 


Examples  of  CPI  Images  of 
Plate-like  Crystals  showing 
Crystal  Edges  (Prism  Face) 


: 


I - 1 

200  urn 

Examples  of  Edges  on  CPI 
Images  of  Irregular  Crystals 


|  | 


Examples  of  Glass  Beads 
Imaged  by  CPI  in  the  Lab 


Fig.  7.  Examples  of  CPI  images  in  SVC  that  show 
edges  (prism  faces)  of  the  larger  plate-like  and 
irregular  crystals,  and  examples  of  images  of  glass 
beads  taken  in  the  laboratory  for  the  sake  of 
comparison 
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Crystal  Habits  for  Particlaa  >  <5  um  [N  >  222) 


Fig,  8.  Comparison  of  crystal  habits  for  (top)  all 
crystals  and  (bottom)  crystals  >  65  pm  in  SVC. 
Particle  habits  were  classified  manually  based  on 
CPI  images. 

The  average  bulk  properties  shown  in  Table  1 
were  derived  by  computing  a  time-weighted 
average  from  each  of  the  19  WB-57F  legs  in  SVC 
(i.e„  the  sum  of  the  number  of  1-Hz  data  points  in 
each  leg  times  the  leg  average  divided  by  total 
number  of  data  points  in  all  legs).  The  average 
particle  concentration  (66  L'1),  extinction 

coefficient  (0.009  km"1),  effective  radius  (8,82  pm) 
and  IWC  (0.055  mg  m"3)  are  very  close  to  values 
derived  from  graphs  shown  in  McFarquhar  et  at 
(2000).  McFarquhar  et  al.  measured  IWC  values 
from  0.001  to  0.1  mg  m"3,  and  show  an  average 
value  of  about  0.05  during  a  descent  through  SVC 
in  their  Fig.  2.  They  measured  a  value  of  about  1 1 
pm  for  effective  radius  on  the  same  flight  Based 
on  cloud  lidar  system  (CLS)  measurements  from 
the  NASA  ER-2  research  aircraft,  McFarquhar  et 
al.  show  a  mean  value  of  i  =  0.045  for  optical 
depth,  based  on  an  average  CLS  cloud  thickness 
measurement  of  0.47  km.  Using  WB-57F 
porpoising  maneuvers  to  estimate  SVC  thickness 
in  CR-AVE  yields  a  range  from  about  0,5  to  1  km, 
which  translates  to  an  optical  depth  of  0.045  to 
0.09 

Jensen  and  Pfister  (2004)  simulated  formation 
of  TTL  cirrus  along  trajectory  curtains  throughout 
the  tropics  using  a  detailed  cloud  model. 
Temperature  profiles  in  the  model  were  taken  from 
meteorological  analyses  with  wave-driven 


perturbations  superimposed.  The  model  tracks  the 
growth  and  sedimentation  of  thousands  of 
individual  ice  crystals  within  each  column  of  air 
(see  Jensen  and  Pfister,  (2004)  for  details).  The 
simulated  mean  ice  concentration  (55  L"1)  and 
effective  radius  (7.3  pm)  agree  remarkably  well 
with  CR-AVE  measurements  (T able  1 ).  Note  that 
if  we  do  not  include  gravity  wave-driven 
temperature  perturbations  in  the  model,  then  the 
simulated  ice  concentrations  (18  L'1)  are  much 
lower  than  the  observed  values  (66  L'1 ). 


Fig.  9.  PSD’s  of  (top)  number  concentration, 
(middle)  area  and  (bottom)  mass  based  on  -1800 
km  of  WB-57F  2D-S  and  CAPS  measurements 
collected  in  SVC.  Green  trace  is  interpolated 
between  5  pm  CAS  and  10  pm  2D-S  size  bins. 
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3.4  Large  particles 


Fig.  10.  Comparison  of  PSD’s  from  drizzle 
observed  in  the  RICO  project  using  (green) 
standard  2D  processing  technique,  (red)  only  in- 
focus  2D-S  images  and  (blue)  images  that  have 
been  resized  based  on  the  Korolev  (2007) 
algorithm, 

Table  1.  Bulk  properties  of  SVC  observed  by  the 
WB-57F  on  4  missions,  and  particle  concentration 
and  from  numerical  simulations  (Jensen  and 
Pfister  2004).  Means  and  standard  deviations  (o) 
are  an  average  of  8738  1-Hz  averaged 
measurements  collected  by  WB-57F  from  19  legs 
in  SVC.  Max  and  Min  are  leg-averaged 
maximums  and  minimums  of  the  19  legs.  _ 


Mean 

a 

Max 

Min 

Particle 
Concentration 
(No.  L'1):  WB-57F 

66.0 

30  8 

188.8 

22.5 

Simulation 

55 

Particle 

Concentration  > 

65  nm  (No.  L*1) 

0.004 

0.017 

0.08 

0.00 

Rrf  (pm):  WB-57F 

8.82 

2.44 

16  7 

5.51 

Simulation 

7.3 

Extinction  (km1) 

0.009 

0.011 

0  063 

0,002 

IWC  (mg  m"3) 

0,055 

0.098 

0.503 

0.012 

The  ice  particles  observed  in  SVC  during  CR-AVE 
extend  to  larger  sizes  than  previously  reported  by 
Heymsfield  (1986),  McFarquhar  et  al.  (2000)  and 
Thomas  et  al.  (2002).  There  were  35 
measurements  out  of  a  total  of  8738  1-Hz  2D-S 
measurements  that  contained  ice  particles  >100 
nm,  and  the  maximum  infocus  particle  size 
observed  is  165  pm.  The  maximum  length  of  the 
largest  CPI  image  is  151  pm.  However,  the 
previous  investigations  were  incapable  of  properly 
measuring  these  larger  particles  (with  the  possible 
exception  of  the  replicator  used  by  Heymsfield 
1986).  McFarquhar  et  al.  do  mention  that 
collateral  studies  in  1973  with  the  Aeromet  Learjet 
near  the  Marshall  Islands  using  PMS  probes 
report  ice  particles  up  to  140  pm  in  thin  cirrus 
between  14  and  15  km  (45,000  and  49,000  ft),  but 
these  clouds  are  much  lower  than  the  SVC  that 
exists  at  the  base  of  the  tropopause. 

While  the  larger  particles  occasionally 
observed  during  CR-AVE  do  not  have  a  significant 
impact  on  bulk  microphysical  properties,  they  do 
raise  the  question  of  how  such  relatively  large 
particles  can  exist  near  the  tropopause.  Jensen  et 
al.  (2007)  use  back  trajectory  analysis  to  show  that 
the  crystals  did  not  originate  in  convection.  The 
large-scale  vertical  ascent  velocity  in  the  TTL  has 
been  estimated  to  be  low,  on  the  order  of  a  cm  s'1 * 
or  less  (Heymsfield  1986).  Growth-sedimentation 
calculations  suggest  the  crystals  cannot  grow  to 
these  sizes  before  falling  out  of  the  supersaturated 
TTL  layer  unless  the  ice  supersaturations  are  on 
the  order  of  100%  or  more  (Jensen  et  al.  2007). 
This  gives  rise  to  discussions  of  discrepancies  in 
results  from  different  water  vapor  measurement 
technologies  that  are  used  in  the  TTL. 

Balloon-borne  measurements  of  frostpoint 
temperature  in  the  TTL  began  in  1949  (Barrett  et 
al  1950)  and  have  been  considered  the  only 
reliable  technology  until  the  recent  development  of 
tunable  diode  and  quantum  cavity  lasers  (May 
1998;  Sayres  2006). 1  The  growth  of  large  ice 

particles  and  how  this  process  is  affected  by 
measurements  of  water  vapor  in  the  TTL  is  the 
subject  of  a  companion  paper  (Jensen  et  al. 
2007).  In  this  paper  we  consider  some 
relationships  between  ice  particle  size  and  water 
vapor  measurements  from  the  ICOS. 


1  A  workshop  that  focused  on  discrepancies  between  balloon- 

borne  frostpoint  and  loser- based  measurements  of  water  vapor 
and  other  aspects  of  the  TTL  was  convened  by  NASA  in 

November  2006 
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Figure  11  shows  a  plot  of  maximum  particle 
size  versus  ICOS  water  vapor  for  8,734  1-Hz 
measurements  of  water  vapor  in  the  TTL.  The 
ICOS  archived  measurements  were  used  only 
when  the  instrument  had  stabilized  after  passing 
through  regions  with  high  water  concentrations, 
and  missing  1-Hz  ICOS  measurements  were 
interpolated  using  existing  measurements.  The 
relationship  between  maximum  particle  size  and 
relative  humidity  with  respect  to  ice  (RH**)  shows 
a  subtle  trend  whereby  the  largest  particles  are 
only  observed  in  regions  at  high  RH„.  For 
example,  all  of  the  SVC  particles  z  120  pm  are 
observed  at  RH**  z  125%.  However,  there  was 
no  observable  trend  when  maximum  particle  size 
was  plotted  against  total  particle  concentration  and 
mean  particle  size. 

Fig.  12  shows  a  time  series  with  RHra, 
temperature,  total  particle  concentration  and  total 
particle  surface  area  on  22  January  2006  where 
appreciable  ice  is  only  found  at  RH„  >  150%,  and 
even  higher  values  of  RH**  correspond  with  the 
lowest  values  of  both  total  particle  concentration 
and  total  particle  surface  area.  Extensive 
measurements  of  water  vapor  in  mid-latitude  cirrus 
at  temperatures  <  -40°C  in  both  the  northern  and 
southern  hemispheres  show  that  more  than  75% 
of  the  measurements  in  both  hemispheres  range 
between  90%  <  RH,„  <  110%  (Ovarlez  et  al. 
2002).  The  authors  suggest  that  the  predominance 
of  RHm,  a  100%  measurements  is  a  result  of  an 
equilibrium  condition  in  these  cirrus  clouds.  In 
contrast,  Figure  13  shows  that  more  than  half  of 
the  CR-AVE  ICOS  RH**  measurements  in  SVC 
exceeded  160%,  which  Koop  et  al.  (2000) 
suggests  is  the  value  at  which  homogeneous 
nudeation  occurs.  Murphy  et  al.  (2006)  and 
Murphy  (2007  -  personal  communication)  report 
that  TTL  aerosols  are  composed  mainly  of  sulfates 
mixed  with  organics.  The  obvious  differences 
between  the  mid-latrtude  cirrus  data  and  our 
observations  are  the  colder  temperatures  and 
unusual  aerosol  chemistry  in  the  upper  TTL. 
However,  there  are  no  obvious  physical 
mechanisms  that  explain  the  observed  SVC 
particle  sizes  and  shapes. 

4  Summary 

During  the  2006  CR-AVE  project  staged  from 
Costa  Rica,  the  NASA  WB-57F  often  observed 
SVC  in  the  temperature  range  from  -75°  to  -85°  C 
within  about  1  km  of  the  tropopause  between 
about  0  and  10°  N  Latitude.  The  average  ice 
particle  number  concentration  (66  L'1),  extinction 


coefficient  (0.009  km'1),  ice  water  content  (0.055 
mg  rn3),  effective  radius  (8.82  pm)  and  optical 
depth  (0.045  to  0.09)  are  similar  to  values 
computed  from  in  situ  measurements  over  the 
Marshall  Islands  by  McFarquhar  et  al.  (2000). 
However,  CPI  images  of  the  ice  particles  observed 
in  CR-AVE  were  markedly  different,  with  84%  of 
the  particles  being  quasi-spherical,  and  particles 
larger  than  about  65  pirn  exhibiting  disk-shaped 
and  hex  plate  characteristics.  In  comparison, 
McFarquhar  et  al.  (2000)  report  a  50%  mixture  of 
columns  and  trigonal  plates  (based  on  Heymsfield 
(1986)  replicator  data  collected  in  1973). 


Composite  of  4  SVC  Missions 


Fig.  11.  Plot  of  maximum  particle  dimension  from 
1-Hz  2D-S  data  in  10  pm  bins  versus  ICOS  RH„ 


Fig.  12.  Example  of  time  series  of  2D-S  particle 
concentration,  temperature  and  ICOS  RH„ 

Also,  the  size  of  CR-AVE  particles  extended  to 
165  pm,  while  McFarquhar  et  al.  (2000)  found  the 
largest  particles  in  TTL  SVC  were  50  pm,  in  both 
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Particle  Concentration  (No  L  \| 


Number  of  1  Hz  Occurrences 


the  1973  WB-57F  study  and  an  investigation  in 
1993  using  a  Learjet. 


measurements  in  SVC  when  RHra  exceeded  100%. 

The  apparent  difference  in  particle  size 
between  the  Heymsfield  (1986)  data  and  recent 
CR-AVE  measurements  may  be  due  to  limitations 
of  the  PMS  1-D  probe  used  in  the  earlier  studies 
(Lawson  et  al.  2006).  However,  the  difference  in 
the  shape  of  the  particles  appears  to  be  real  and  is 
not  readily  explained.  One  possible  difference  is 
that  very  high  relative  humidity  (in  excess  of 
200%)  observed  by  the  WB-57F  during  CR-AVE 
was  not  present  in  the  1973  case  reported  by 
Heymsfield  (1986).  Based  on  a  numerical  model, 
Jensen  et  al.  (2007)  show  that,  for  crystals  with 
aspect  ratios  of  6:1  or  less,  3  to  A  ppmv  water 
vapor  (i.e. ,  about  200%  to  250%  RH,„  at  -85°C)  is 
necessary  to  grow  ice  particles  with  sizes  on  the 
order  of  100  urn  within  about  500  m  of  the 
tropopause.  Also,  chemistry  measurements  near 
Costa  Rica  in  2004  and  during  CR-AVE  suggest 
that  the  TTL  is  composed  mainly  of  sulfates  mixed 
with  organics  (Murphy  et  al.  2006;  Murphy  2007  - 
personal  communication).  Nucleation  and  growth 
of  ice  particles  at  -85  °C  in  the  presence  of  sulfates 
mixed  with  organics  and  very  high  relative 
humidity  are  mechanisms  that  need  to  be 
investigated  more  thoroughly,  perhaps  with  the  aid 
of  cloud  chamber  experiments. 
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Appendix  A  Calculation  of  particle  aspect  ratio 
from  CPI  images 

Two-dimensional  particle  imagery  (e.g.,  CPI 
images)  limits  the  amount  of  three-dimensional 
information  retrievable  from  any  given  image. 
However,  the  aspect  ratio  (an  essentially  two- 
dimensional  quantity)  can  be  estimated  from  some 
CPI  images  that  view  the  crystal  from  a  fortuitous 
angle.  Estimates  of  particle  aspect  ratio  (i.e.,  ratio 
of  dimensions  of  basal  to  prism  faces)  based  on  a 
CPI  image  from  the  CR-AVE  data  set  that  reveals 
three-dimensional  qualities  are  discussed  in  this 
Appendix. 

Figure  At  shows  a  magnified  and  rotated  view 
of  the  crystal  third  from  the  left  in  the  top  row  of 
images  in  Fig,  7.  The  crystal  appears  to  be  a 
hexagonal  plate  viewed  from  such  an  angle  that  a 
basal  face  and  at  least  two  prism  faces  are  seen. 
Assuming  hexagonal  symmetry  an  estimate  of  the 
view  direction  and  thus  of  the  aspect  ratio  can  be 
made.  Drawn  over  the  image  in  Fig.  A1  are  yellow 
lines  approximating  the  locations  of  the  six  edges 
of  the  basal  face.  The  relative  (i.e.  arbitrary  units) 
lengths  of  these  lines  are  indicated  in  yellow  as 
well.  Opposite  edges  of  the  basal  face  are  similar 
in  length  as  expected  by  symmetry. 

Let  A  represent  the  average  length  of  the 
shortest  pair  of  basal  edges  (i.e.  39.5),  fi 
represent  the  average  of  the  mid-length  pair  (i.e. 
86.5),  and  C  the  average  of  the  longest  pair  of 
edges  (i.e.  99.5).  Let  <J>  represent  the  angle 
between  the  normal  to  the  crystal’s  basal  face  (i.e. 
the  c-axis)  and  the  view  direction  and  9  the  angle 
between  the  red  line  in  the  plane  of  the  basal  face 
and  the  projection  of  the  viewing  direction  onto  the 
same  plane.  Let  L  represent  the  true  length  of  the 
edges  (in  the  same  arbitrary  units).  Basic 
trigonometry  assuming  a  hexagonal  plate  yields: 

A1/!}  =  eos2(4>)cos2(©)  +  sin:(0) 

fi 2//,2  =  oos2  (4>)  cos2  (60  -  ©)  +  sin 2  (60  -  0) 

C3/i2  =  cos2  (<P)cos2  (60  +  0)  + sin  2(6O  +  0) 


Solving  these  equations  for  the  values  of  A,  B,  and 
C  specified,  yields:  L  *  106,  6  *  8°,  and  O  *  70° 


Fig  A1.  Exploded  view  of  CPI  image  of  ice 
particle  shown  in  Fig.  7.  See  text  for  explanation 
of  lines  and  numbers. 


The  pink  line  represents  the  thickness  of  the 
plate  from  the  view  angle.  It’s  length  (36),  in  the 
same  arbitrary  units,  must  be  divided  by  sin(&)  to 
obtain  an  estimate  of  the  actual  thickness  (TH). 
The  aspect  ratio  is  then  the  ratio  of  2 L  I  TH  or 
about  5.5:1.  The  largest  uncertainty  is  in  the 
choice  of  the  thickness  (pink  line).  Figure  7  shows 
that  the  CPI  optics  can  blur  edges  making  them 
look  thicker  than  they  are.  If  this  is  causing  an 
exaggeration  of  the  thickness  then  we  may  be 
underestimating  the  aspect  ratio.  Therefore,  we 
conclude  that  the  aspect  ratio  of  this  crystal  is  on 
the  order  of  6:1. 
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Appendix  B  Uncertainty  in  the  size  of  ice 
particle  images 
B,1  2D-S  Images 

For  2D-S  images  that  appear  to  be  in  focus  (Le., 
all  pixels  are  black,  there  are  no  white  pixels  in  the 
center),  there  are  three  sources  of  uncertainty. 
Each  are  considered  in  turn  below: 

L  Uncertainty  due  to  the  effective  pixel  size  of  the 
photo  diodes  in  the  sample  volume.  For  the 
2D-S  the  effective  pixel  size  is  10  pm. 
it,  Uncertainty  in  the  calibration  of  effective  pixel 
size  due  to  uncertainty  in  the  laboratory 
calibration  of  glass  beads  used  to  determine 
effective  pixel  size.  Based  on  discussion  with 
the  bead  manufacturer  Duke  Scientific,  this 
uncertainty  is  about  5%  of  the  image  size. 

Hi  An  image  can  be  somewhat  out  of  focus,  thus 
enlarged,  and  still  appear  in-focus  (all  image 
pixels  are  black;  no  white  pixels  in  the  center). 
Using  the  tables  found  in  Korolev  (2007),  the 
result  depends  on  the  size  of  the  image  rather 
than  being  a  fixed  value,  or  fixed  percentage  as 
discussed  above  in  items  i  and  //  Estimates  of 
the  maximum  error  are  shown  as  function  of 
image  size  in  Table  B1,  based  on  the  following 
discussion. 

Table  B1  shows  calculations  of  actual  image 
sizes  (D0)  for  six  apparent  image  sizes  (Dm  =  50, 
75,  100,  125,  150,  175  pm)  that  are  slightly  out  of 
focus,  but  not  sufficiently  to  produce  a  white 
(Poisson)  spot  in  the  middle  of  the  image.  For 
item  ill,  the  error  is  always  in  the  direction  of  an 
overestimation  of  the  size.  Therefore  we  find  the 
maximum  overestimation,  minimum  size,  possible 
by  first  assuming  the  effect  of  (1)  above  has 


caused  a  10  pm  overestimate  of  the  apparent  size, 
and  then  assuming  it  is  also  out  of  focus  as  much 
as  possible  and  still  not  have  a  spot.  The  3rd 
column  in  Table  1  is  the  ratio  of  the  spot  size  j(Dt)] 
to  apparent  size,  (from  Korolev  2007)  Using  D%  of 
15  pm  yields  a  (normalized  distance  from  the 
focal  plane  from  Korolev  2007  Table  B1  column  1) 
and  the  ratio  of  apparent  size  to  actual  size  (from 
Korolev  2007  Table  B1  column  2  ,our  column  5), 
from  which  D0  is  computed.  The  maximum 
underestimation,  or  maximum  possible  size,  is 
estimated  as  Dm+10. 

B.2  CPI  images 

Unlike  the  2D-S  probe,  the  CPI  is  designed  to 
maximize  detection  and  imaging  of  only  infocus 
images.  It  is,  however,  possible  for  CPI  images  to 
be  out  of  focus,  but  unlike  the  2D-S  probe  that  has 
only  one  gray  level  (i>e.,  black),  the  CPI  images 
particles  with  256  gray  levels  This  significant 
increase  in  shading  allows  the  eye  to  readily 
identify  CPI  images  that  are  infocus,  especially 
when  the  images  have  relatively  sharp  corners 
and  regular  structures,  such  as  the  plates  shown 
in  Fig  7  in  the  text.  Thus,  uncertainty  in  the  size  of 
infocus  CPI  images  reduces  to  the  pixel  size  (2.3 
jim)  plus  the  amount  of  blur  that  results  from 
optical  aberrations.  Based  on  laboratory 
evaluation,  the  amount  of  blur  is  estimated  to  be 
equivalent  to  about  one  pixel,  so  the  total 
maximum  uncertainty  in  CPI  images  that  are  in 
focus  is  estimated  to  be  about  5  pm. 


Table  B1:  Dm  is  the  observed  size  along  the  array,  D0  is  the  smallest  possible  image  size 
(i.e.,  the  actual  particle  when  in  focus)  and  Dm  +10  Is  the  largest  possible  image  size 
given  errors  discussed  in  the  text. 


Dm 

Dm  -  10 

D^(Dm-10) 

Z<j 

(Dm-10)/Do 

D0 

Dm+10 

50.000 

40.000 

0.375 

2.420 

1.166 

34.299 

60  000 

75.000 

65.000 

0.231 

1.810 

1412 

46.030 

85000 

100.000 

90.000 

0.167 

1.058 

1142 

78.800 

110.000 

125.000 

115.000 

0.130 

0.900 

1.238 

92.694 

135.000 

150.000 

140.000 

0.107 

0.709 

1191 

117  570 

160  000 

175.000 

165,000 

0.091 

0.580 

1143 

144403 

185.000 
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